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1 Introduction 

The description of the long-time, low temperature, collective dynamics of CQudensed 
matter systems is one of the central aims of quantum many body theory .BBB Nu- 
merous successful theories have been developed for a variety of materials, and we 
can broadly separate them into two (not entirely distinct) categories: 
Order parameter dynamics: Often the system is in or near a phase with a bro- 
ken symmetry. This can be used to define an order parameter which, along with 
globally conserved quantities, obey equations of motion which define the collective, 
long-time dynamics. The simplest example of this was the van Hove theory, which 
evolved into a sophisticated description of the hishly non-trivial dynamics near fi- 
nite temperature second order phase transitionso 

Quasiparticle dynamics: The elementary excitations above a particular ground 
states are identified, and a transport equation for their collision and interactions 
provides a description of the low frequency response functions. Landau's Fermi 
liquid theoryu is the familiar example of such a theory. 

One common feature of both classes of models is that the dynamic equations con- 
trolling the lowest frequency dynamics are ultimately classical. They do, never- 
theless, provide a description of quantum systems at low temperature. Quantum 
mechanics does play a fundamental role in determining the nature of the ground 
state, its excitations, and of effective coupling constants in the equations of motion, 
but the equations are finally expressed in terms of classical, collective degrees of 
freedom. For the order parameter dynamics, the justification for a classical de- 
scription is that the typical relaxation frequency is smaller than ksT /h (where T is 
the absolute temperature), while in the quasiparticle models we require a collision 
frequency smaller than ksT /fi. 
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In this paper, we will review recent studies of finite temperature dynamics near 
second-order quantum phase transitions We will find that the two classes of classi- 
cal models discussed above do apply over a significant portion of the phase diagram. 
However, we will also discuss the novel quantum critical region,El 'Vilifire both de- 
scriptions are inadequate. The fundamental property of this regionuu is that the 
characteristic frequency is a number of order unity (which is usually universal) 
times kBT/h: so neither classical description is evidently applicable. A theory of 
quantum-critical dynamics can, of course, be obtained from an exact solution of the 
critical field theory describing the phase transition: we will see an example of this in 
Section ^ and the results yield considerable insight. However, such exact solutions 
are rare, and cannot be relied upon for a general theory. Vaciaus perturbative ap- 
proaches to quantum critical dynamics have been developed,Bll3 and these usually 
rely upon effective classical models of order-parameter or quasiparticle dynamics as 
a point of departure. 

We will begin our discussion in Section ^ by considering the quantum Isiug 
chain. An essentially exact description of the long time dynamics is availablalil 
in all the different low temperature regions of the phase diagram of this model, 
and this will allow us to describe and distinguish the essential characteristics of 
the order-parameter, quasiparticle, and quantum-critical dynamics. Section |^ will 
consider a quantum critical point in two dimensions between two superconductors 
with distinct pairing symmetries. The critical theory involves interacting fermionic 
and bosonic excitations, and we will discuss approximate theories for the quantum 
critical dynamics and their possible relevance to recent photo-emission experiments 
on the cuprate superconductors. Finally, in Section ^ we will discuss the quantum 
critical dynamics in simplified models of correlated electrons in the presence of 
quenched disorder, and their relationship to experiments on the heavy fermion 
superconductors. i_. 

The reader is also referred to another recent reviewllj which covers additional 
topics on quantum phase transitions presented in my talk at MB- 11. 

2 Quantum Ising chain 

We consider one of the simplest models which displays a second-order quantum 
phase transitions; its phase diagram is nevertheless rich enough to display all the 
dynamic regimes we wish to distinguish. The Hamiltonian of the quantum Ising 
chain is 

j \l>Q / 

where (t|, are Pauli matrices on a chain of sites j, Je (> 0) are short-ranged 
exchange constants, and g (> 0) is the transverse field. The operator JJj aj com- 
mutes with iJ/p-and generates its global Z2 Ising symmetry. We will consider the 
phase diagramtll of Hi as a function of g and T, shown in Fig |l|. The global Z2 
symmetry is spontaneously broken at T = for g < gc, where gc is a critical value 
of order the Ji. There is a second order phase transition at g — g^ T — 0, and the 
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ure 1. Phase diagram of the quantum Ising chain H[ in (111). There is long range order, with 
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3), only for g < Qc and T = on the shaded line. The dashed lines represent cross overs 
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at A ~ ksT. Quasi clas sica l models for the long times dynamics in the lower left (Section 
and low er r ight (Section 2.2) regions are available. The quantum critical dynamics is discussed in 
Section 2.3. 



symmetry is restored elsewhere in the phase diagram. We describe the different 
regimes of Fig ^ in the following subsections. 

2. 1 Order parameter dynamics 

First consider Hj at g ~ Q. Here Hj is the classical Ising model and all states are 
known exactly, as is the T > partition function. The ground states are doubly 
degenerate, I T)i or Y\j I and have long-range ferromagnetic order in o-|. 
All excited states can also be written down exactly, and the lowest energy ones are 
domain walls between the two ordered states 

...|T>j-2|T>,-iIT>jUUiU>j+2U>,+3..., (2) 

whose energy is A = 2 X]^>o ^"^f- a-bovc the ground state. The entire set of excited 
states can be built out of multiple domain wall states. At a small T > 0, there is 
an exponentially small density p ^ g-A/fc^T thermally excited domain walls. A 
simple calculation shows that the presence of these excitations is sufficient to destroy 
the long range order beyond an exponentially large correlation length ^ — l/2p. 

Now move to a small g > 0. All the states are now perturbed by quantum 
fluctuations, but the qualitative picture remains the same at long length scales. 
There is still ferromagnetic order in the ground state, and a broken Z2 symmetry, 
with the long range correlation 

lim {a]at)T=o = Nl (3) 

\j-k\^(X} 

but A^o is no longer unity as it is in the classical model; instead quantum fluctuations 
reduce Nq with increasing g until it ultimately vanishes at g = gc- This subsection 
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will restrict attention to g < gc- The degenerate manifold of the domain wall 
states (H) now broaden into a band of single particle states with dispersion Sp = 
A+AgsvcL'{p/2)+0{g'^) as function of momentum p; we assume unit lattice spacing) 
and have now redefined A to be the true energy gap including all corrections. The 
resulting motion and collision of the thermally excited domain walls is responsible 
for the long time reiaxationai dynamics of the ferromagnetic ordeiiparameter. This 
picture was used to develop a simple solvable quasiclassical theoryllil for the dynamic 
correlations of the order parameter in the g < gc, low T region of Fig |l|, and the 
results are most conveniently expressed in terms of the dynamic structure factor 
S{p,uj): 

(4) 

At T = 0, the long-range order (|3|) leads to the elastic Bragg peak 

S{p,uj) = {2^fNl5{p)5{uj) + ... (5) 

where the ellipsis represents contributions at frequences cj > 2A associated with 
the creation of two or more domain walls by the external probe. At T > the 
quasiclassical order parameter dynamics broadens the delta functions in and 
leads to 

S{p,uj) = / dx die-^(P^-"*)i?(a;, t) + . . . (6) 



where R{x,t) is a reiaxationai function. The form of R{x,t) for general x, t is 
involved, but its basic features can be understood by looking at two simple limits; 

i?(a:,0) = e-l^l/« ; ^ = 2 f ^e'^'"'-'^ 
i J-TT 27r 

So R decays exponentially on a characteristic spatial scale ^, and a characteristic 
temporal scale, the coherence time r^; the delta functions in (||) therefore broaden 
at r > to a momentum width of order 1/^, and a frequency width of order 
1/T(^. The expression ^ contains an exact result foiLx^ at low T: the quantum 
Ising chain, and the closely related dilute Bose gas,E3Ej are the only interacting 
many body quantum systems for which exact results for such a relaxation rate 
are available. Remarkably, involves only the energy gap A and fundamental 
constants of nature; it also satisfies r^p ^ h/ksT as is required for the applicability 
of the quasiclassical theory. 



2.2 Quasiparticle dynamics 

Next consider the opposite limit of very large g. At g = cx) all eigenstates can again 
be written down simply. There is a unique, paramagnetic ground state Y[j I 
where | — >)j = (| t)i + I l)j)/V2 is an eigenstate of crj. In strong contrast to the 
states at g — 0, equal time correlations of f7| are now non-zero only on-site, and so 
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the spin correlation length is effectively zero. The lowest excited states are obtained 
by flipping a single spin: 

... I ^>,-2| ^>,-l| H. I H, + l| H,+2| H,+3 . . . , (8) 

where | — (| t)j — I i)j)/V2 is the other eigenstate of crj, and these states 
have energy 2g above the ground state. We can now examine corrections to the 
limit case in powers of \/g. The ground state remains paramagnetic but the spin 
correlation length increases as g is decreased. The flipped spin states again 
broaden into a quasiparticle band with Sp = A + 4 X]£>o "^^ sin^(p^/2) + 0{g'^) and 
A = 2g—2 J2i>Q Ji+^id'^)- Keep in mind, however, that the physical interpretation 
of these quasiparticles is entirely distinct from those of Section they are not 
domain walls, but localized spin flips. A crucial property of the quasiparticles is 
that they have an infinite lifetime for a finite range of momenta around p = 0: this 
is simply because energy and momentum conservation prohibit their decay into 
any other states. Consequently, although equal time correlations of cr| decay very 
rapidly in space, there are long-range correlations in spacetime; this is evident from 
the expression for the dynamic structure factor 

S{p, iu) = ttA6{uj - Ep/n) + ... ; T = 0, (9) 

associated with quasiparticle pole, as infinite range temporal and spatial correla- 
tions are needed for the Fourier transform to conspire to have a pole. Here the ellip- 
sis represents multiparticle contributions at frequencies u > 3A. The quasi-particle 
residue A is finite in the paramagnetic phase and decreases as g is decreased. 

At T > 0, there will again be an exponentially small density of these thermally 
excited quasiparticles. A quasiclassical transport equation for the collisions of these 
quasiparticles can be written down and, remarkably, solved exactly. This solution 
shows that the quasiparticle collisions broaden into a Lorentzian 

^(P'^) = 7 it(2ln/ ^2 - (10) 

The expression for the quasiparticle-width, turns out to be identical to that 

in (0), where A is now the energy gap to the flipped-spin quasiparticles. Note that 
the characteristic collision time is again larger than Ti/ksT ^ as is needed to justify 
a classical dynamical model. 



2.3 Quantum critical dynamics 



The quasiclassical models in Sections 2.1 and provide an essentially exact de 



scription of the long time dynamics as T — > for any fixed g, other than at the 
critical point g = gc- As illustrated in Fig |^, this failure broadens into a wide 
quantum critical region at larger T, which is, in principle, easily accessible in ex- 
periments. 

We begin our discussion of spin correlations in this region by first considering 
the single point g = g^ T = 0. The energy gap vanisheSphere as A ^ \g — gd- 
Also as g /' gc, the ferromagnetic order parameter vanishesEj as A'o ~ {gc — gY^^, 
and so it is not sensible anymore to thinkJji terms of domain walls. Similarly, for 
g \ ffc, the quasiparticle residue vanisheaid as A ^ (g — gc)^^^, and so the flipped 
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spin quasiparticles are also not well defined &i g = Qc- A less intuitive description 
of the long time dynamics must be developed here. The fundamental property of 
the critical point that enables a description of the dynamics is that of scale and 
conformal invariance. One signal of this invariance is the long distance behavior of 
the equal time correlation 

K'^^) = ^Z^ 5 -9 = 3- ^ = 0' (11) 

where Z is a non-zero, non-universal number. Remarkably, the powerful technology 
of conformal invariance allows one to reconstruct the exact long distance and long 
time spin correlation function at T > using as input only the form (^) (such 
'wizardry' is unfortunately not possible in higher dimensions). The answer is more 
easily expressed in terms of the dynamic susceptibility, x(p, w), which is related to 
S{p,uj) in (0) by the usual fluctuation dissipation theorem: 



r(7/8)c'^/-^ \ Z \16 AwkBT/hJ \16 AirkBT/h 

27/47r3/4r(l/8) ) {kBT/ny/'i^ fW . u + cp \ . u-cp 



16 AnkBT/hJ \16 AnkBT/h 

where c is the non- universal velocity of excitations at the critical point. This is a 
rather cumbersome expression whose functional form is not terribly intuitive. Its 
structure is nevertheless rather simple, and this becomes clear from the follow- 
ing expression which provides an excellent fit to (^2|) over a wide window of low 
frequencies or wavevectors (this is also the region over which the spectral density 
Imx(p, cj)/w has significant weight): 

x(o,o) 

i-iiu;/rH)+P^e-i^/^i) 

The relaxation rate Tn is given by 



X(P,i^) = -, ,,,.,r._^ , 2f2 7-7— T2- (13) 



TT \ kflT 



LJ = 



.(2tan-)^ (14) 

where we used ( p^ ) in the first equation. As claimed in the introduction, this 
characteristic frequency is a universal number times kBT/h. We determined uji 
and ^ by fitting ( p^ to ( p^ over the range < u!,cp < 2kBT/h (the fit was 
essentially perfect) and obtained the best fit values 

LOi = 0.795kBT/h ; f = 1.280hc/kBT. (15) 



As expected, all time and length scales are determined by T alone. Notice that ( p_3D 
has the structure of th&.response function of the conventional van Hove dynamics 
of an overdamped field.EJ What is novel here the characteristic damping frequency 
is not an adjustable phenomenological parameter depending upon the coupling to 
some heat bath, but a universal rate determined by the absolute temperature and 
fundamental constants of nature. 
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At large frequencies, h\uj\ S> fc^T, and wavevectors, fic\p\ 3> fcsT, the expression 
2|) reduces to 

, , r(7/8)27/Vc3/4 Z 

^^^'") = f(T78^ (cV-^^)^/« ^"^^ 

This represents the spectrum at the gapless critical point g = gc,T = Q] as expected, 
the complex frequency plane has no quasiparticle pole, but a branch-cut originating 
at Lu = ±cp associated with the continuum of critical excitations. For practical 
purposes, it is worth noting that the spectral weight in this critical region is rather 
small, and experimental observations at T > will be dominated by the relaxation 
spectrum in (p^. 

3 Transition between BCS superconductors 

This section reviews quantum critical dynamics in a two-dimensional model which 
mayj-be of relevance to the cuprate high temperature superconductors. Observa- 
tionta of quantum critical behavior in -the spectral function of the fermionic nodal 
quasiparticles stimulated our studytj'tZl of a wide class of possible models which 
could yield a suitable quantum critical point. It turned out to be rather difhcult 
to satisfy the needed requirements: a quantum critical point below its upper crit- 
ical dimension, with nodal quasiparticles a primary degree of freedom. Only two 
such candidates were found, both involving transitions with changes in the pair- 
ing symmetry of BCS superconductors in two dimensions: the transitions from a 
dx^_y2 superconductor to either a (1^2 _y2 -f idxy or a (1^2 _y2 + is superconductor. 
We considered the first of these more likely on microscopic grounds and will review 
its theory here; the theory for the second is closely related. We also note that our 
interpretation of the spectra of the cuprates does not require that such a quantum 
phase transition be actually present on the physical axis as a function of doping; 
all that is required is that a quantum critical point be nearby in the generalized 
parameter space. 

Recent tunnelling measurementalj have provided support for a transition from 
a (1^2 _y2 superconductor at optimal doping to a 6,^2 _y2 + idxy superconductor in the 
overdoped regime, at least in thin films of Yi_yCayBa2Cu307_2;. This sequence of 
transitions may be understood on the basis of some early theoretical work. While 
there is strong theoretical evidence for dj.2_j,2 superconductivity in lightly doped 
antifecmmagnets,!!^ studies of the Hubbard model in limit of very large hole den- 
sityl23'Ell (low electron density) showed instead an instabilitv to d^y superconduc- 
tivity. Interpolating between these limits, we can expectll3 that a dj.2_j^2 -|- idxy 
superconductor appears an intermediate phase. 

As we are interested only in universal critical properties, we can develop a theory 
using the simplest phenomenological model which displays the required transition. 
So we consider the Hamiltonian 

Htj = ^ ekcl^Cka + Ji ^ Si • Sj + J2 ^ Si - Sj, (17) 

A; (ij) nnn ij 

where Cja- is the annihilation operator for an electron on the square lattice site 
j with spin a =T)ij ^ka is its Fourier transform to momentum space, is the 
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dispersion of the electrons, 



Sja = ^c],0,c,.. (18) 



with cr" the Pauh matrices, and Ji, J2 are first and second neighbor antiferro- 
magnetic exchange interactions which induce the superconductivity. We apply the 
standard BCS theory to Htj: this will yield an adequate description of the low 
temperature properties, except in the vicinity of the quantum critical point. The 
BCS Hamiltonian is 

Hbcs = ^fc^^cfc- ^A.(sV]+/ia - sVI+a.t) + 

-yE'^^(4T^Ui-^]i4+p,T) + h-, (19) 

where the first summation over ^ is along the nearest neighbor directions x and 
y, while the second summation over v is along the diagonal neighbors x + y and 
—a; + y. To obtain d^2_y2 and d^y pairing, we choose Aa; = —Ay = A.j.2_y2, and 
Ax+y — —A^x+y — Axy We summarize our choices for the spatial structure of the 
pairing amplitudes (which determine the Cooper pair wavefunction) in Fig ^. The 
values of Aj.2_j,2 and A^y are to be determined by minimizing the ground state 
energy: 



EbCS = J\\Ax2-y2\ +J2|A2;j' 



el + \Ak\^-ek 



(20) 



47r2 

where the pairing amplitude in momentum space is 

Ak = JiAx2-y2 (cos kx — cos ky) + 2JiAxy sin kx sin ky. (21) 
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Notice that the energy depends upon the relative phase of ^^i^yi and IS.xy'- this 
phase is therefore an observable property of the ground state. 

The minimization of Ebcs was carried out in Rcf. ^ For small J2/J1 the 
ground state was a dx^-yi superconductor with A^2_j^2 ^ and l^xy — 0. Above a 
critical value of J2/ Ji there was a continuous transition to a state with /!\x2_y2 ^ 0, 
A^j^ ^ 0, and 

arg(A^j,) = arg(A^2_j^2) ± 7r/2. (22) 

This is a (1^2 _y2 +idxy superconductor. The choice of the sign in ( p2|) is associated 
with the breaking of time reversal symmetry, and this suggests that the transition 
can be described by a Z2 Ising order parameter. As usual, we can associate the 
Ising order with a real, scalar field (f) (the coarse-grained value of the local Ising 
order) which we can identify here by = iA^y, in the gauge where Ax2_y2 is real. 
However, the critical theory of the transition is not simply that of the quantum 
Ising model in 2+1 dimensions i.e. not that of the two-dimensional version of 
the model considered in Section |^. The fermionic quasi-particles are also primary 
critical degrees of freedom, and the required quantum field theory couples the Ising 
field (p to fermionic fields representing the lowest energy Bogoliubov quasiparticles. 
We will not explicitly write down this field theory here, and instead refer the reader 
to another recent reviewli3 by the author. 

The finite temperature phase diagram of Hf j in the vicinity of the critical point 
is sketched in Fig ^. We discuss below the nature of the fermion Green's function 
in the different regimes below the superconducting critical temperature T^. 

We can write the fermion Green's function in the superconductor compactly 
using the Nambu notation. We define the spinor field ^'^ = (cfc|,cLj,|). Then the 
retarded '^k Green's function G{k,uj) can be written as 

hG-\k,uj) = ?iw-efcT^ -f Re(Afc)T^ -Im(Afc)T^ - I](fc,cj) (23) 

where r^'^'^ are Pauli matrices in Nambu particle-hole space. The above expression 
follows from the BCS theory (|l9|), with the self energy S representing the effects of 
interaction between the quasiparticles and any additional collective modes. 

Consider first the low T region with J2 < J2C in Fig ||. Here the ground state 
is well described by the BCS theory of the (ij.2_j,2-wave superconductor. We can 
compute the effects of collisions between the fermionic quasiparticles in very much 
the same spirit as in Landau's Fermi liquid theoryfl the main difference is that 
instead of a whole Fermi surface, we now have isolated 'nodal' Fermi points at 
four points in the Brillouin zone determined by the solution of ek = 0, Afc = 0. 
Consequently, the density of states for low-energy scattering of quasiparticles is 
even smaller than in a Fermi liquid, and the quasiparticles remain well-defined 
excitations. Right at the Fermi level, and at low T, these collisions lead to E ~ 
iT^, which is smaller than the result ~ iT"^ in a conventional Fermi liquid. A 
quasiclassical description of the fermionic quasiparticles therefore remains valid, in 
the same framework used in Landau's Fermi liquid theory. 

Next, we move to the opposite low T region in Fig ^ with J2 > J2C and dx2-y2 + 
idxy superconductivity. Here time-reversal symmetry is broken with (</>) ^ 0. The 
theory of the quasiparticles is very similar to that in the small J2 region above 
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Figure 3. Phase diagram of Htj ■ As in Fig hi, dashed lines represent crossovers, but the full lines 
are true thermodynamic phase transitions. There is a quantum critical point at T = 0, J2 = J2C1 
but now it extends into a line of second order phase transitions at terrmeratures Tj > 0. The 
critical theory for the T = quantum critical point was reviewed in Ref. ha and involves an Ising 
field and fermionic degrees of freedom. The Tj line is in the universality class of the classical 
two-dimensional Ising model. The ground state is superconducting everywhere below Tc, with 
<i^2_y2 + idxy superconductivity in the shaded region, and d^2_„2 superconductivity elsewhere. 



The crossoveir^mndaries and Tj approach the T : 
best estimatecj'Lj for the exponent zu is 1.00. 



: 0, J2 = J2c critical point linearly because the 



with one very important difference. The equations Sk — 0, Afe = now have no 
soiution for any k and so the nodal points have been gapped out. There are still 
low-energy fermionic excitations near the position of the original nodal points, but 
these appear only above an energy gap A. As in Section |[ this energy gap mnishes 
as we approach the critical point as A ( J2 — J2cY'^ , and the best estimatecj'O for 
the critical exponent is zv « 1.00. As long as T ^ at any fixed A > 0, the density 
of thermally excited quasiparticles is exponentially small, and so S ^ e^^l"^ . So 
again a quasiclassical transport theory of the quasiparticle collisions applies. 

Finally, we turn to the intermediate temperature quantum critical region in 
Fig m in the vicinity of the J2 ~ J2c critical point. Here the physics is rather 
similar to the corresponding region in Section 2.3, but we no longer have the benefit 
of an exact solution for the spectral function. Strong fluctuations of the Ising order 
parameter 0, whose condensate gapped the nodal points in the J2 > J2C region, now 
lead to strongly damped fermionic excitations. Loosely speaking, rapid temporal 
oscillations between dx2_y2 -\- idxy and dj.2_y2 — idxy pairing lead to broad spectral 
functions. A theory for the low frequency form of the fermion spectral function 
was developed in Ref. |l6|, using an expansion in (3 — d) (where d is the spatial 
dimensionality) and phenomenolgical ansatzes similar to (^3|) for the frequency and 
wavevector dependence of the spectral function of the critical excitations. At low 
frequencies near the Fermi level, the result can be written in the form (compare 
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@ and (|l|)) 

G~\k, u) = Z-iT"/ {hu - EkT' + AfcT" + ihTn) , (24) 

where 77/ is a critical exponent, Z is a non-universal, non-singular prefactor, and we 
can take real because (0) = in the quantum critical regiDn. The estimates for 
the exponent 77/ are 77/ « {3 — d)/llAii the (3 — d) expansionjlj and rjf ~ l/{3n'^N), 
with = 2, in the 1 /N expansionB3 Th e da mping rate Tr has the same remai|kable 
universal structure we found in Section [2T^ . Just as in (p^, we now flncO'tJ that 
Tr is universally related to the absolute temperature, with 

= 0.581%^ (25) 

in the (3 — d) expansion. 

It is important to keep in mind that (^J) is a low frequency form and does not 
hold at frequencies which arc much larger than ksT/h. In the latter regime we 
obtain a result which is similar in structure to that obtained by taking the large 
u! limit of the exact result (|lj) of the quantum Ising chain. Near one of the nodal 
points we find 

G{p,.) = ZCj (26) 

where p measures momentum deviation from one of the nodal points and the p 
co-ordinates have been rotated by 45 degrees from the k co-ordinates, C/ is a 
computable universal number, and c is a velocity. In principle the velocities appear 
before the and should bje-jdifferent, but they become equal in the asymptotic 
region near the critical point .E3 The spectrum ( |2^ ) is the analog of (p^), and its 
most important property is, of course, that it does not have quasi-particle pole, but 
only a branch cut representing the continuum of critical excitations. 

4 Disordered spin systems 

The plethora of intermetallic "heavy fermion" compounds have provided a fertile 
ground for the study for the study of quantum phase transitions and the associ- 
ated quantum critical region at finite temperature. Useful reviews of the current 
experimental jaiid theoretical status have been provided recently by Coleman and 
collaborators .E3 Most of these systems are near a zero temperature magnetic order- 
ing transition of some type, and this has allowed investigation of quantum criticality 
using a variety of probes. Among the most thoroughly investigated compounds is 
CeCue-ajAua;: neutron scattering and thermodynamic measurements^ show con- 
vincing signs of universal critical behavior, with a characteristic frequency scale of 
order ksT /fi and scaling of dynamic, response functions as a function of hjjj /ksT . 
This behavior is not compatibles^ Eil with spin-density wave theories,c3^EZl which 
consider perturbative corrections from interaction between collective paramagnon 
modes, and the search for the appropriate quantum-critical model is an active topic 
of current research. An important open question in these investigations is the im- 
portance of quenched disorder to the observed critical behavior. Theories for the 
onset of spin glass order in metallic systems have been proposed (see Ref. pSl for 
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a review), and glassy dynamics has been observedS in UGu^-i^Pdj;, a compound 
which also displays huj/ksT scaling in neutron scattering.L^ More unexpectedbi, 
the glassy behavior persists at compositions which are nominally stoichiometricEil 

Here we shall review a class of theories which have their origin in studies of sys- 
tems with stmng quenched disorder. Analyses, of disordered models of frustrated 
spin systemsH'Ej and doped Mott insulatoraij have shown that a useful starting 
point for the critical theory is a seemingly simple model which we shall describe 
below. It consists a single quantum spin interacting with a 'bath' of low energy 
collective spin excitations which represent a mean-field description of its environ- 
ment. Generalized models have recently been argued—to apply to systems in two 
dimensions even in the absence of quenched disorder 

We consider the following single spin quantum partition function, which defines 
a well-posed, quantum mechanical problem in its own right: 



l(T))dT 



Z = j Vn{T)5{n^ {t) - I) cyip (^-iS Ar{n{ 

+ dT dr'DiT - T')n(T) • n(T')^ 



(27) 



Here r is imaginary time which runs periodically from to /3 = h/kBT, n(r) is 
the orientation of the quantum spin with nir + 13)= n(T), and S its spin quantum 
number with 2S = integer. The first term in Z is the Berry phase of the spin, 
with (n(r))dr equal to the area of the spherical triangle on the unit sphere with 
vertices at n(T), n(r -I- c?t), and a fixed (but arbitrary) reference direction Hq. The 
function D{t) (also periodic with period j3) is the retarded self-interaction of the 
spin and represents the infiuence of the spin cmiironmcnt. The model Z was intro- 
duced and examined in jB-laxge N expansion,L3 and useful rcnormalization group 
analyses were presented. E3' til In its application as a mean-field theory of bulk lat- 
tice models, the solution of Z usually has to be supplemented by a self-consistency 
condition for D{t): the simplest of these is D{t — t') = (J/?i)^(n(r) • n{T'))z, 
where J is an energy scale for the exchange interacliniis-Jput more complicated 
self-consistency conditions have also been considered.oE3'E3 

In a quantum critical system, we may expect power-law correlations of all ob- 
servables at T = 0. So it is of particular interest to examine Z for the case where 

Dir) - (28) 

for large |t| at T = 0. Non-tciyial spin correlations emerge for a > 0, and the 
rcnormalization group analysif£3EZI proceeds in powers of a. Although this is an 
interacting quantum field theatipjid it is not possible to compute the flow equations 
exactly, it was recently shownHEj that the structure of the quantum theory allowed 
one to deduce some results to all orders in a; in particular it was shown that for 
the case where D{t) obeyed (|2F 



(n(r).n(r'))^^^-^. (29) 
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Figure 4. Plots of the dynamic local susceptibility ImxL('^), in (|3l[), and the local dynamic 
structure factor (measured in neutron scattering experiments) Sl{llj) = 2ImxL — e"'''^''*^-^) 
at a = 1/2 and some fixed temperature T. Note that Imxi('^) is an odd function of ui, while 
5i(a;) is not an even function. At large |ti;|, ImxLi'^) ~ sgn(ii;)/|ti;|^~". 



This is a non-trivial result,ilSi3iEii and depends crucially on the the Berry 
phase term in Z. Comparison of (^9|) and ( p8| ) with the particular self-consistency 
condition mentioned above selects the value [O. |^ J— jand this is the value that has 
appeared frequently in physical applications.tHlESBEJ 

The result ( |29| ) allows us to directly compute the observable local dynamic spin 
susceptibility, Xi('-^) in the quantum critical region; we define this quantity by the 
Fourier transform in imaginary frequency 

XLiiuJn) = f dTe'^-^n{T) ■ n{T'))z. (30) 
Jo 

While the result ( ]29| ) is exact at T = for the model obeying (p^), the form of the 
exact result at T > is not kno-p|i|_.However, a solution for the T > dynamics 
is possible in the large TV limit the right hand side of (|2^) is replaced by 
[ttT/ sin(7rr|r — t'|)]", and its Fourier transform yields the result 



Im[xLH] -r"-isinhf^ 



\2kBT 



a . huj 



2 27rkBT 



(31) 



This constitutes the complete description of the relaxational quantum critical spin 
dynamics: note again that the function depends universally on Tilj /ksT . At the im- 
portant value a = 1, ( |3l| ) reduces to the simple result Im [xL(t^)] ~ tanh(/iaj/2fcBT). 
We show a plot of the spectral function (|3l]) at a = 1/2 in Fig 0. 
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The result ( pl| ) is also that expected if the local spin dynamics obeys a confor- 
mally invariant field theory in 1+1 dimensions .e3 There is no reason for this to be 
generally the case for Z, and (|3l|) has been shown to hold only in a large N limit; 
we expect corrections to the scaling fmiction of w/T at higher orders, and comput- 
ing these remains an important challonga. IJpmparisons of (^) with experimental 
data on UCus-a^Pda; have been made;Ej'CJ'c2l while some higher frequency features 
are reasonably captured, the agreement is poor faL.?iw <ti ksT, and it has been 
proposed that this is due to the effects of disorder .e3 



5 Conclusions 



We have presented here a series of paradigms of dynamic response functions at 
finite temperatures near a quantum critical point. The behavior discussed here is 
characteristic of critical points below their upper critical dimension. As the number 
and precision of experimental measurements of such critical points increases, we 
hope that the accuracy of theoretical predictions will be eventually be sufficient to 
permit a quantitative confrontation between theory and experiment. 

We have not discussed here the dynamics of systems above their upper-critical 
dimension: in these cases the dynamics usually remains quasiclassical even near the 
critical point, and a van-Hove like theory with a relaxation rate dependent upon 
microscopic details provides an adequate description. 
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